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ABSTRACT 

In this paper we show how the rescattering of CMB photons after cosmic reionisation can give 
a significant linear contribution to the temperature-matter cross-correlation measurements. 
These anisotropies, which arise via a late time Doppler effect, are on scales much larger 
than the typical scale of non-linear effects at reionisation; they can contribute to degree scale 
cross-correlations and could affect the interpretation of similar correlations resulting from the 
integrated Sachs-Wolfe effect. While expected to be small at low redshifts, these correlations 
can be large given a probe of the density at high redshift, and so could be a useful probe of 
the cosmic reionisation history. 

Key words: cosmic microwave background — large-scale structure of universe. 



1 INTRODUCTION 

Power spectrum measurements of the CMB temperature and the 
large scale distribution of matter have long been important tools 
for probing cosmology, but only recently has the study of correla- 
tions between these important probes been proved to be as useful. 
Indeed, the absence of strong cross-correlations was instrumental 
in showing tha t the CMB anisotro pies we observe were of cosmo- 
logical origin dBennett et al.ll 19931) . However, the relatively small 
cross-correlations that do exist give us an important way of con- 
straining those subdominant temperature anisotropies that are cre- 
ated relatively locally. 

Two of the most important local sources of temperature 
anisotropies contributing t o the cross-correlation s are the integrated 
Sachs-Wolfe (ISW) effect dSachs & Wolfell967t) and the Sunyaev- 
Zeldovidich (SZ) effect dSunvaev & ZeldovictJ [l97(ih . Both of 
these effects provide means of following the growth rate of struc- 
ture at late times, and through it can tell us about the recent acceler- 
ated expansion of the Universe. The ISW fluctuations are generated 
when the gravitational potential begins to evolve at late times due to 
the influence of dark energy; this is a linear ef fect and the expected 
cross -correlations are on large angular scales dCrittenden & Turokl 
On the other hand, SZ fluctuations arise from scattering 
from hot ionised electrons in clusters, and so lead to smaller an- 
gular cross-correlations. 

Recently, evidence for both of these effects has been 
seen in cross-correlation stud ies between the CMB anisotropies 
seen by the WMAP satellite dBennett et ai]|2003l ; lllinshaw et all 
2006) and various surveys of large scale structure. Large an- 
gular correlations consistent with the ISW effect have been 
observed in cross-correlations with radio, infrared, x-ray and 
optical data, particularly that from the Sloan Digital Sky 
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Survey (SPSS) dBoughn & Crittenden! I200I; iNolta etai] 120031: 
' If 



l Afshordi etai] 120031: IFosalba & Gaztafiagaj 120031; IScranton et al 
| 2003 |: IFosalba et alj|2003l; IPadmanabhan et all 120041 ; ICabre et al. 
booalRassat et alj|2006l: iGiannantonio et alJl200fj) ~The SZ effect 
is usually observed in targeted cluster observations, but evidence 
for it has also been seen in cross-correlations studi es, for example 
betwe en WMAP and the 2MASS infrared survey dAfshordi et al.l 
120031) . 

While the ISW and SZ effects appear to be the dominant 
sources of cross-correlatio ns (once more mundane foregrounds 
have been excluded, e.g. Gia nnantonio et all d200fj|) ). there are 
other local sources of correlations which are potentially impor- 
tant to understand in order to accurately interpret the observa- 
tions. On large scales, it has recently been shown that grav- 
itational magnification can project local inhomogeneities onto 
higher redshift surveys, causing the matter density at those red- 
shifts to appear more correlate d with the ISW anis otropies than 
would be expected otherwise dLo Verde et al. 2006). On smaller 
scales, effects l ike the kinetic SZ, the clo sely related Ostriker- 
Vishniac effect dOstriker & Vishniad[T986T) and the Rees-Sciama 
effect dRees & Sciamal 19681) could tell us much about the evolution 
of structure , and particularly help probe it s velocity on these scales 
dlliev et alj 120061; iMcOuinn et all 12005k ISchaefer & Bartelmannl 
l2005l;IStebbinsll2006h . 



In this paper, we examine another possible source of cross- 
correlations. Like the kinetic SZ and Ostriker- Vishniac effects, 
it results from Doppler scattering off of moving electrons, lead- 
ing to CMB anisotropies with the same frequency dependence 
as the primordial anisotropies. Here we focus on the Doppler 
anisotropies resulting from the electron velocities arising at linear 
order, which can appear on large angular scales; these could po- 
tentially affect the interpretation of the ISW effect. The impact of 
this effect on the CMB temperature power sp ectrum is well un- 
derstood, where it is known t o be subdominant {b odelson & Jubasl 
1 19931 : ISugivama etai] 1 19931 : iHu et alj [19941 ICoorav & Hulll999l) . 
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lAlvarez et al ] ( 120051) showed that reionisation can produce a sig- 
nificant correlation between the 21cm HI radiation and the CMB; 
here we focus on its impact on CMB -galaxy cross-correlation mea- 
surements, where the effect can be comparable to the ISW at high 
redshifts, and if unaccounted for, would bias the estimation of pa- 
rameters. 

In section 2, we discuss reionisation and outline the late time 
linear contributions to the CMB auto- and cross-correlation spec- 
tra; in sectio n 3 we present the results coming from a version of 
CMBFAST dSeliak & Zaldarriagall 19961) . an extension of modifi- 
cations made in lcorasaniti et alj d2005t) . We discuss the prospects 
for observing the effect in section 4, before drawing conclusions in 
section 5. 



2 THE EFFECTS OF REIONISATION 



2.1 Reionisation history 



Cosmic reionisation is currently thought to be caused by the UV ra- 
diation emitted by the first luminous objects, and is experimentally 
constrained by the optical depth to electron sc attering of CMB pho - 
tons, given by WMAP as K r = 0.092 ± 0.030 dSpergel et ai]|2006t) . 
It is also constrained by measurements of the Gunn-Peterson 
troughs in the Lyman-a part of the spectra of distant quasars from 
the SDSS dFan et alj[200ll : IWhite et alj|2003h . These suggest that 
the inter-galactic medium (I GM) was fully ionis ed out to a red- 
shift £ = 6.10 ± 0.1 5 dGnedin & Far] 120061) . As noted by 
IShull & Venkatesanl d200?t) . the scattering up to this redshift ac- 
counts for nearly half of the total observed optical depth. 

Precisely how reionisation happened prior to this is a topic of 
some debate; the minimal assumption is that the universe became 
completely ionised very quickly at a single redshift (around a red- 
shift of z r = 11 to be consistent with the optical depth constraint.) 
However, the process could have been more complex. In the fol- 
lowing, we focus on two other models: a "double step" model with 
ionisation fraction brought first to 1/2 at zi = 15 and then to 1 at 
Z2 = 6, and a parametrisation of the double reionisation model by 
iSr] (2003), which again has two distinct phases at z\ = 15 and 
Zi = 6, but with an evolving ionisation fraction between them. Fig- 
ure Q] shows the visibility function g(z) as a function of redshift in 
these different reionisation scenarios. Here, the visibility function 
gives the probability that a photon last scattered at a given redshift, 
and is related to the optical depth by g(r) = n e cr r ae~* (z) . 



2.2 CMB anisotropics from reionisation 



When reionisation is introduced, a second peak appears in the vis- 
ibility function, corresponding to the restored coupling between 
photons and matter at late times. Because the CMB photons can 
now again scatter off free electrons, their properties will thus be 
altered by the temperature, potential, and velocity of the scatterer, 
exactly in the same way it happened at the last scattering surface 
at z = 1100; these however are suppressed by the small frac- 
tion of photons that are rescattered at this time. At linear order , 
there are three main effects of reionisation dSugivama etaill 19931 ; 
lHu&Dodelsonll200lh : 

• Damping of anisotropics on small scales, 

• Secondary anisotropy production after reionisation, 

• Additional polarisation on large scales. 
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Figure 1. Visibility function g(z) for the best fit WMAP third year ACDM 
model, with and without reionisation. Double reioni sation models are also 
shown, for a double step case and for the|c3 J2003I) scenario. 



In addition, other nonlinear effects such as the Ostriker-Vishniac 
effect and patchy reionisation will be relevant at smaller angular 
scales ({ > 1000). The small scale damping is by far the most 
dominant effect, affecting the temperature power spectrum at small 
scales (£ » 10) by a factor ~ e~" . The polarisation production is 
also well studied and is responsible for breaking the degeneracy 
between the optical depth and th e spectral index o f the primordial 
perturbations in the WMAP data dPage et alj|2006l) . 

Here we focus on the impact of the anisotropies generated af- 
ter the photons rescatter, which are largely dominated by the mo- 
tion of the scattering electrons. Along any given line of sight, one 
expects electron velocities to be moving towards us at some red- 
shifts and moving away from us at other redshifts, leading to a can- 
cellation of the Doppler effects. However, the scattering probabil- 
ity is not uniform, causing the Doppler effect to be dominated by 
the redshifts soon after reionisation. For this reason, there is a net 
anisotropy produced in the CMB temperature, which is correlated 
with the additional polarisation produced that this time. The greater 
the gradient of the scattering probability, the less the cancellation 
and the larger the anisotropies. 

Figure [2] shows the total anisotropy power spectrum with- 
out and with reionisation and the single late time contributions in 
the reionisation case. The damping on small scales is the biggest 
contribution, and we can see the small role played by ISW and 
velocities, while the combined effect of density and gravitationa l 
pot ential perturbations is neg ligible, as found by iHu et alj d 19941) 
and iDodelson & Jubasl dl993l) . In a scenario without reionisation, 
gravitational effects like the ISW would be the only source of 
anisotropies at late times. The Doppler contribution, while always 
subdominant, is actually larger than the ISW contribution on suffi- 
ciently small scales. At even smaller scales, the thermal SZ effect 
is important, although it can be distinguished from the other effects 
taking advantage of its frequency dependence. 

More interesting is to consider how the Doppler contribution 
might be correlated with the density. Consider matter falling into 
a large over-density at some redshift zo', the matter on the far side 
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Figure 2. Total temperature anisotropies power spectrum for the best fit 
WMAP third year ACDM model, with (thick solid line) and without (thin 
solid line) reionisation. The Doppler anisotropies produced after reionisa- 
tion (dashed line) are significantly smaller than the ISW anisotropies (dot- 
dashed line) on very large scales, but can be comparable on degree scales. 
At even smaller scales, the therm al SZ effect becomes relevant (the dot- 
ted theoretical curve is obtained by Schaefer et al. (2004) using the W-body 
simulation method bv lSpringel & Hernquistl 1200 II) ). although it can be dis- 
tinguished from the other effects taking advantage of its frequency depen- 
dence. The reionisation model is the default single step. 



(z > Zo) will be travelling towards us, while the matter on the near 
side (z < Zo) will be travelling away from us. Scattering from both 
sides will contribute to temperature anisotropies, but because the 
scattering is more likely at higher redshifts, the electrons moving 
towards us will be more likely to scatter. Thus the over-density will 
be associated with a temperature hot-spot, and in the same way 
under-densities will be associated with CMB cold-spots. 

In the following, we will show that the effect of this secondary 
rescattering, though always negligible for the temperature power 
spectrum, is present and can be important in the temperature- 
matter correlations at high redshift, and in particular is comparable 
to the magnification bias effect; therefore this effect must be taken 
into account to produce precise cross-correlation predictions. 



Figure 3. Total temperature-density correlation power spectrum for the best 
fit WMAP third year ACDM model and a matter visibility function centred 
in z = 3, with reionisation. The reionisation model is the single step and the 
galaxy bias is set to 1 . 



gauge, where the numerics are usually performed, their peculiar 
velocities will be very small. 

We usually work in the harmonic space using a Fourier expan- 
sion for quantities like the velocity 



v(x,t) = ^kv(k,r)e' kx , 



(2) 



(where we have assumed the velocity is irrotational.) Following 
Ma & Bertschingeil dl995l) , we can relate the Newtonian gauge ve- 
locities to those in the synchronous gauge by, 



v = v, - ika, 



(3) 



where a = , h and rj parameterise the metric degrees of freedom 
in synchronous gauge and dots refer to derivatives with respect to 
the conformal time. Since we know the synchronous gauge veloc- 
ities will be small (v s m 0), the Doppler term will be dominated 
by the term proportional to a. We can relate this to the dark matter 
density fluctuation using its conservation equation: 



2.3 Power spectra 

We next describe how to calculate the power spectra for these 
secondary Doppler temperature anisotropies and their cross- 
correlations with the density. These will be proportional to the frac- 
tion of photons that are scattered and to the velocity of the scatter- 
ing electrons v relative to the observers line of sight n: 



A r (n) 



-r 



dr g(r) v(x, t) • n, 



(1) 



where x = (tq— T)n, To is the present conformal time and t, is taken 
to be the time just prior to the beginning of reionisation. While the 
physical effect will be gauge invariant, this expression is implicitly 
in the Newtonian gauge. By the epoch of reionisation, the electrons 
will have fallen into the dark matter wells and will be effectively 
at rest with respect to the dark matter particles; in the synchronous 



6 = — = -k?a + 3n - -lea. 
2 ' 



(4) 



The Einstein equations give i] cc v,; so, assuming we can neglect 
the synchronous velocities, we find: 



(5) 



Going back to the anisotropies, and defining as usual fi = k ■ n, 
we have 



A r (n) = - fVjy 

Jt: i. 



S(k,r) 



(6) 



which can be integrated by parts; dropping the surface terms where 
the visibility function is small, we find 



A r (n) = - jH dr J] e**^ I L [g(T)<S(k, t)J . 



(7) 
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Figure 4. Temperature-density correlation power spectrum for the best fit WMAP third year ACDM model in function of z, with the set of galaxy selection 
functions defined in Eq. (TS) (redshift tomography). The Doppler effect is compared with the ISW. We see the different evolution of the three different 
reionisation histories described in Figure [T] to every change in the sign of the visibility function derivative corresponds a change in the sign of the cross- 
correlation due to the Doppler effect. The galaxy bias is always set to 1 . 



If we assume linear theory for the growth of density perturbation 
so that 6(k, t) = D(r)S(k, 0), we can then expand the exponential 
in terms of spherical harmonics to show 



A r (n) = 



- f ° dr J] <5(k, 0)^ ^ tMKjo - T))K ; „,(n)y ( ;(k) 

Jr i k K lm 



X -[g(T)D(T)\ 



(8) 



From this expression, we find the harmonic coefficients for the 
Doppler anisotropies to be 



aj m = Y<5(k,0)47n'y; m (k)A' r (fc), 



where 



Aj(k) 



1 n 

Jr, 



dTj t (k(r -T)) — [g(T)D(T)\ 



(9) 



(10) 



The expectation of the square of these coefficients gives the 
Doppler power spectrum 



c ' T =U 



dkk 2 P(k)\Aj\ 2 , 



(11) 



where P(k) is the matter power spectrum today. 

Following similar arguments, we can find simi lar equations 
for th e anisotropies in the matter density (see e.g. iBoughn et al.l 
dl997l) ). Specifically, we find 



<, = y^(k,o)Wy; m (k)A; ? w 



but instead 



A*(k)= drb g (r)D{T)W{r) Mk(r -t,), 



(12) 



(13) 



where b g (r) is the possibly evolving galactic bias and W(j) the mat- 
ter selection function normalised to unity. We can then use these 



to derive the matter auto-correlation and the matter-temperature 
power spectra 



Cf = - f P(k)k 2 \Al(k)\ 2 dk 
Cj g = ^ f P(k)k 2 A] {k)A s ,(k)dk. 



(14) 



For the following discussion, we implement these in a numerical 
Boltzmann code, a modified version of CMBFAST. 

At sufficiently small angular scales, it is possible to simplify 
the calculations using the small angle approximation. Here, we can 
integrate t he complete e volution Eq. dl It in an approximate analyt- 
ical form jLimbej[l954h : 



Cj T 



c 



Tg 



2 Cdrt l+Ml Ud 
n J r 2 \ r j\dr 
2 C dr 



g(r)D(r) 



a + i/2) 2 

d ~ 



(15) 



dr P (l±ltll bgW{r)D(r)dr 



g{r)D{r) 



a+i/2) 2 



where r = t -t is the conformal distance. However, on large scales 
this approximation will not be valid, and therefore throughout we 
compute all our results using the full-sky Boltzmann code. 

We can see from Eq. l !15b how this effect depends on the cos- 
mological parameters: the derivative of the visibility function g will 
bring a dependency on the reionisation history, as we will see in 
more detail in section [3j on the other hand, the derivative of the 
growth factor depends on the matter parameters Q.,„, tr 8 , and also on 
the dark energy equation of state. However, in typical models the 
visibility function will change more dramatically than the growth 
function, and will dominate the effect. 

Here we have focused on the cross-correlation with the mat- 
ter density; the cross-correlation with the 21 cm line depends on 
instead o n the neutral hydro gen density which can be more com- 
plicated ^Alvarez et al.ll2005l) . If reionisation were uniform, peaks 
in the matter density would be associated with a higher neutral hy- 
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drogen density; however, ionisation is expected to be produced by 
the early structures formed in the highest peaks, causing the peaks 
in the matter density to have less neutral hydrogen. Depending 
on the competition between these effects, the CMB-21 cm cross- 
correlation can have the opposite sign from what is derived above 
dAlvarez et ai]|2005l) . However for the discussion below, we will 
always assume the tracers are positively correlated with the total 
matter density. 



2.4 Double scattering 

As noted bv lCoorav &"h3 dl999f) . double scattering processes can 
produce higher order corrections to the Doppler signal. These cor- 
rections have a power spectrum which, for scales smaller than the 
width of the visibility function, can be approximated by 



_l r 



dTj,[k(r - T)]g(T)k(T)D(T). 



(16) 



Their effect is small, always < 5% at the peak, and we will therefore 
neglect it, this being anyway a conservative choice. 



3 RESULTS 

3.1 Low redshift signal 

The expected cross-correlations will depend on many factors, but 
the biggest are the depth of the survey and the assumed history of 
reionisation. To begin, we look at the amplitude of the effect at 
low redshifts. Figure [3] shows the expected cross-correlation with 
a broad selection function at mean redshift 2=3. This assumes a 
selection function of the form, 



3z 2 4 
W,(z)=--e \ 



(17) 



where z = 1.41z , which is a good approximation to the distribu- 
tion of objects in a flux limited survey, without any imposed cut 
in redshift. For comparison, we show the ISW cross-correlation for 
the same survey depth. For this depth, the effects are actually com- 
parable, though this is in part because the ISW effect peaks around 
z — 1 and is already dying off at these redshifts. 

3.2 Reionisation history dependence 

The low redshift signal is largely independent of assumptions of 
the reionisation history, since we believe the Universe to be totally 
ionised at that time. However, the signal will increase as we probe 
higher and higher redshifts and will provide us with a way to probe 
the visibility function evolution. To study the time evolution of this 
effect, we use narrower selection functions corresponding to the 
imposition of redshift cuts o n a survey; to be spec ific, we use a 
redshift slicing proposed by iHu & Scrantonl d2004) for an LSST 
type survey. There they assumed the total selection function could 
be divided into many narrower bins at different redshifts. The z'-fh 
selection function was assumed to be 



Wi(z) oc W,(z) x \ erfc 



(i-l)A- 



r(z)V2 



- erfc 



(A- 



r(z)V2 



(18) 



where A = 0.8, <t(z) = 0.02(1 + z) and every slice is normalised to 
unity. To study the evolution out to the time of reionisation, we will 
consider i = 1, 19, corresponding to mean redshifts z, from 0.5 
to 15. 



We can see how the effect evolves with redshift in Figure [4] 
The ISW effect is dominant at low redshifts (where dark energy is 
dynamically important) and becomes negligible at high redshifts. 
The Doppler contribution largely depends on the time derivative 
of the visibility function. At late times, once the Universe is fully 
reionised, the visibility function decreases smoothly as the uni- 
verse expands and the density of scatterers decreases; the result- 
ing cross-correlations evolve relatively slowly. However, the pic- 
ture can change dramatically once we reach redshifts where the 
ionisation fraction changes, which can drastically alter the visibility 
function. This is particularly the case if we had a survey sensitive 
to the beginning of the ionisation history, which for the WMAP 3 
best fit model is at z — 11. There, the sign of the correlation be- 
comes negative and the amplitude can be large if ionisation occurs 
quickly. 

This dependency on the visibility function makes it a sensitive 
probe of the reionisation history. In Figure [4] we compare the sig- 
nal for different histories: the single step reionisation, a simple dou- 
ble step rei onisation, and a more complicated reionisation based on 
|C^(2003). At low redshifts, all models are fully ionised and have 
equivalent cross-correlations. When the histories diverge (z > 4), 
the cross-correlation can be reduced or change signs twice if we 
are able to see two epochs of the visibility function increasing. 

The most prominent feature is actually the negative cross- 
correlation which occurs at the primary epoch of reionisation. As- 
suming a single step process, the amplitude of this effect depends 
on two factors: the total optical depth and how long it takes to 
ionise the universe. The total optical depth determines the redshift 
of reionisation as well as height of the peak of the visibility func- 
tion. However, since the time derivative of the visibility function 
comes into the integral, the duration of reionisation is also impor- 
tant. 

The full impact of these factors is most easily seen in the tem- 
perature spectrum, which sums the anisotropies from all redshifts. 
Increasing the optical depth raises the signal and shifts the spectrum 
to slightly smaller scales. On the other hand, shortening the dura- 
tion of reionisation allows the resolution of velocities on smaller 
scales, and so adds small scale power to the anisotropies. 

3.3 Comparison with magnification 

It is interesting to compare the late Doppl er effect with the grav - 
itational magnification effect described by iLo Verde et al.1 {2006); 
while no additional temperature fluctuations are created by the 
magnification, this effect can make the low redshift density dis- 
tribution appear to be in a higher redshift sample, thereby causing 
it to be more correlated with the ISW anisotropies than one might 
expect. The amplitude of these will depend on how the number den- 
sity of objects depends on the flux cutoff. Since no new anisotropies 
are created, the effect introduces covariances between the low and 
high redshift cross-correla tions . 

ILo Verde et al.l d2006l) showed that ignoring the magnification 
effect can bias the value of the dark energy equation of state which 
would be inferred from the cross-correlation measurements. We 
show in Figure [5] that the additional cross-correlations from the 
Doppler effect are of a similar magnitude to the gravitational mag- 
nification effect, so one would expect a similar kind of bias in the 
equation of state if they were ignored. But unlike the magnifica- 
tion effect, the Doppler correlations represent true new correlations 
at these redshifts, so the covariance between the higher and lower 
redshift measurements will be much smaller. We studied the bias 
that would be obtained by ignoring the Doppler correlation (but 
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Figure 5. Comparison of the magnification and Doppler cross-correlations. 
We show the ISW effect alone (blue, dot-dashed), the effect of adding 
corrections from cosmic magnification (red, dashed) and all three ef- 
fects together (black , solid). Here we have reproduced the calculations of 
lLoVerde et ailfeOOd) using the same cosmology and assumptions about the 
galaxy samples (mock catalogue I) and dN/dz. 



Figure 7. The total signal to noise for ISW and rescattering Doppler effect, 
for the single step reionisation history. The different lines represents differ- 
ent redshift cuts, respectively from top to bottom: no cut and cuts at redshift 
10, 8, 6. The insert shows the total signal to noise (without redshift cuts) for 
the velocity effect. At high {, cross-correlations from non-linear effects will 
become important. 
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Figure 6. The inferred value of the dark energy equation of state pa- 
rameter wobtained ignoring the effect of reionisation, for the samples of 
lLoVerde et alj )2006l) . The input model had w = -1. The top panel shows 
the result for the single redshift bins, while in the bottom panel we present 
the effect of adding up to the i-th bin. 



in cluding magnification ) with a Fisher matrix analysis, as done 
bv lLoVerdeetalJj2006h . The inferred values for w obtained with 
marginalising over the other parameters Q. m , Q.t, n s , h and erg with 
priors errors of 5% are shown in Figure [6] The overall inferred w 
differs from the "true" one that we put into the calculation by al- 
most 2a. 



4 OBSERVABILITY 

The ISW effect itself is difficult to detect from cross-correlation 
measurements, principally because of the presence of large pri- 
mordial CMB anisotropies originating at a redshift of z ~ 1000. 
These fundamentally limit the significance which the the cross- 
correlations can be observed to a signal to noise level of 7 - 10. 
It is worth exploring if the late Doppler anisotropies are similarly 
limited. 

We first consider the most optimistic picture possible, in which 
we had a full sky map out to a given redshift and were able to 
reconstruct a map of the predicted late Dopple r anisotropies. In this 
case the expected signal to noise is given by dCrittenden & Turokl 
1 1995b 



(2/ + D- 



C 



T-Dop 



(19) 



By looking at the auto-correlations in Figure|2] we see that for large 
(, the Doppler anisotropies exceed the ISW anisotropies. Thus we 
might hope that since the signal to noise weights the large multi- 
poles more, the total signal to noise might be higher than the ISW 
case. This is indeed the case, as demonstrated by Figure [7] When 
we include the whole signal, which includes the cross-correlations 
from the epoch of reionisation, the signal to noise can be much 
larger than for the ISW. For a typical model (k = 0.092) we find 
a total optimal signal to noise of order S /N ~ 20. This can rise to 
as high as 30 with a higher optical depth and shorter duration of 
reionisation, or drop as low as 10 if we take these parameters to the 
other extreme. However, the figure also shows that the much of the 
small angle anisotropies arise at very high redshift, and it would 
take a very deep survey before the significance could exceed more 
than a few. 

The above calculations are for the most optimistic case. A 
more realistic estimate can be found using the calculated cross- 
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Figure 8. Evolution of signal to noise ratio in function of the redshift bins 
for the velocity effect in different reionisation scenarios compared with the 
ISW. 



still be important in the interpretation of the integrated Sachs- 
Wolfe signal. For example, we have shown that they are compa- 
rable to the effect of gravitational magnification. A high redshift 
ISW signal has also been shown to be a potential means of discrim- 
inating dark e nergy models from those in which the laws of gravity 
arc modified dSong et al.| [2006); to test such models correctly, in- 
cluding the cross-correlations from rescattering would be essential. 

Detecting these cross-correlations in their own right will 
clearly be a challenge since it requires deep probes of the matter 
density. One possibility is to use galaxy or qu asars found with a n 
optical survey such as th e DES (for quasars) jAbbot et alj|2005h . 
LSST Jlvson et al.l l2006). and possibly the panoramic survey of 
the SNAP, which is still in the proposal stage, which could probe 
the Universe to a redshift of z = 3. Alter natively, radio instruments 
such a s LOFAR teottgering et alJ|2006l) or the SKA dBlake et al.l 
|2004 iTorres-Rodriguez & Cresj [20071) could potentially go as 



deep or deeper by finding 21 cm HI emission from early galaxies; 
or, if the near-infrared back ground is dominated by the first genera- 
tion o f structure formation ( Santos et al. 2002; Sa lvaterra & Ferraral 
2003), it might be useful as a proxy for the de nsity at high redshifts. 
Finally, as discussed by I Alvarez etaU J2005l) . the 21 -cm emission 
from the neutral gas at reionisation could be observed directly by 
an experiment like SKA, providing a probe of the density precisely 
where the ionisation is changing most dramatically. 



correlations functions for the galaxy selection functions described 
above. In this case, the signal to noise is given by 

In (zi) \ = V + Vfa^r+[cp"p - (20) 

We show this signal to noise ratio for the Doppler contribution com- 
pared to the naive ISW effect in Figure U for the WMAP best fit 
model of quasi instantaneous reionisation and the series of redshift 
tomography proceeding from the matter visibility functions defined 
in Eq. J 1 8t >. We can see how a small but non negligible S /N is pro- 
duced at high redshift, while a much larger signal is generated at 
the deeper reionisation epoch. 

A small but potentially measurable signal to noise is present at 
redshifts z < 7, where we know that we can find galaxies and other 
collapsed density tracers; the signal can be enhanced if the ionisa- 
tion is still evolving at these low redshifts. A much larger signal 
is expected at higher redshifts, but there the best density probe i s 
probably the 21 -cm radiation as suggested bv lAlvarez et al.l f2005h . 



5 CONCLUSIONS 

In this paper, we have focused on the linear contributions to cross- 
correlations arising from the rescattering of CMB photons. These 
arise from scales of order k ~ 0.01/? Mpc~ and are seen on de- 
gree scales. However, it should be emphasised that these are by 
no means the only such sources of cross-correlations; non-linear 
effects such as patchy reionisation, the Ostriker-Vishniac effect 
and high redshift SZ sources all could potentially contribute to 

cross-correlations from the reionisation epoch (see the recent pa- 

I 1 1 

per by Slosar at al. ( 2007)). These cross-correlations will typically 

be at much smaller angular scales. While the linear effect is sen- 
sitive to the bulk properties of reionisation, such non-linear cross- 
correlations could potentially provide useful information on the de- 
tailed physics of the reionisation process itself. 

At low redshifts, these cross-correlations are small but could 
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